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We present a survey of the phenomenological adjustment of the parameters of the Skyrme-Hartree- 
Fock (SHF) model for a self-consistent description of nuclear structure and low-energy excitations. 
A large sample of reliable input data from nuclear bulk properties (energy, radii, surface thickness) 
is selected guided by the criterion that ground-state correlations should remain small. Least squares 
fitting techniques are used to determine the SHF parameters which accommodate best the given 
input data. The question of the predictive value of the adjustment is scrutinized by performing 
systematic variations with respect to chosen nuclear matter properties (incompressibility, effective 
mass, symmetry energy, and sum-rule enhancement factor). We find that the ground-state data, 
although representing a large sample, leave a broad range of choices, i.e. a broad range of nuclear 
matter properties. Information from giant resonances is added to pin down more precisely the 
open features. We then apply the set of newly adjusted parameterizations to several more detailed 
observables such as neutron skin, isotope shifts, and super-heavy elements. The techniques of least- 
squares fitting provide safe estimates for the uncertainties of such extrapolations. The systematic 
variation of forces allows to disentangle the various influences on a given observable and to estimate 
the predictive value of the SHF model. The results depend very much on the observable under 
consideration. 

PACS numbers: 21.10.Dr, 21.10.Ft, 21.10.Gv, 21.60.Jz, 24.30.Cz 



I. INTRODUCTION 

The production of exotic nuclei far off the valley of 
stability is making steady progress at various labora- 
tories around the world and is providing an increasing 
amount of new experimental data on basic nuclear prop- 
erties, see, e.g., [J This is a great challenge for nu- 
clear structure theory at all levels of refinement, from the 
more phenomenological microscopic-macroscopic meth- 
ods (see, e.g., 0) through self-consistent mean-field 
(SCMF) methods (see, e.g., Q) or large shell-model cal- 
culations (see, e.g., Q) up to several ab initio techniques 
employing a given nucleon-nucleon interaction (see, e.g., 
00, Hi)- Further development of all these methods is a 
highly topical task for which a large network of activities 
has been launched recently 0. The present paper aims 
to contribute to the development of SCMF. 

There are several approaches to SCMF from which the 
most widely used are the Skyrme-Hartree-Fock (SHF) 
model [3l , the Gogny force [T^j , and the relativistic mean- 
field model (RMF) jii|,[i3. The genuine nucleon-nucleon 
interaction does not allow an immediate mean-field treat- 
ment. Thus all SCMF models employ effective interac- 
tions which are arranged to provide reliable nuclear struc- 
ture properties and low-energy excitations at the level of 
a mean-field description. That approach has much in 
common with the density functional theory widely used 
in the physics of electronic systems [l^l- The difference 
is, however, that electronic correlations are well under 
control and that reliable electronic energy-density func- 
tionals can be derived from well controlled ab initio cal- 
culations. The nuclear case is much more involved be- 



cause the nucleon as such is a composite particle and 
a nucleon-nucleon interaction is already an approximate 
concept. Thus, nuclear many-body theories have not yet 
reached sufficient descriptive power to serve as direct in- 
put for deducing effective energy-density functionals for 
SCMF. Although there are promising attempts for an ab 
initio derivation the general strategy for construct- 
ing energy functionals for high-quality calculations is to 
deduce the formal structure from principle considerations 
im and to rely on a phenomenological adjustment of the 
model parameters. The present paper aims at a criti- 
cal and thorough survey of the SHF model and the phe- 
nomenological adjustment of its parameters. There is a 
long history of SHF development and optimization (for 
recent reviews see and still need for further intense 

consideration, as can be seen from the lively discussion 
in the literature, recent re-fitting attempts [13, [iHl, ^'^d 
the huge combined effort in 0. This study will point 
to several incompatibilities with the present SHF ansatz 
when trying to cover more than just the basic ground- 
state properties. We are thus not yet at the stage to 
advertise one preferred parameterization. The strategy 
is rather to supply a toolbox of forces produced under 
systematically varying conditions. 

After a brief review of the basic SHF functional and the 
strategy of least-squares fitting, our considerations start 
with an inspection of the input data for the adjustment 
of the model parameters. The idea is to select those 
ground-state observables (energies, radii, surface thick- 
nesses) which are expected to be well described by a pure 
mean-field model. The limits to a mean-field description 
are set by the ground-state correlations (GSC) stemming 
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from low- lying collective quadrupole modes, for recent 
smveys see [l^, HO, l2l[. Taking the general trends of 
GSC as worked out in [2l|, we determine a set of obscrv- 
ables which have only small correlation effects and which 
we use as basis for the fits of mean-field models. We then 
will investigate in detail how predictive such fits can be. 
To that end we consider series of fits with additional fea- 
tures added as constraint. These features arc quantified 
in terms of well known nuclear matter properties (NMP) 
like incompressibility if, effective mass m*/m, symme- 
try energy agym, and sum-rule enhancement factor k 
We study systematic variations of these NMP. Further 
information from giant resonances (GR) is invoked to 
define the optimal values for the NMP more precisely. 
Finally, we apply the set of SHF parametcrizations thus 
obtained to a study of the predictive value for more de- 
tailed properties like, e.g., neutron skin, isotope shifts of 
charge radii, and super-heavy elements. 



Pp and pn- In detail they read 



II. THE SKYRME ENERGY FUNCTIONAL 



The goal of a nuclear mean-field theory is to describe 
the many-body system exclusively in terms of a set of 
single-particle wavefunctions together with the BCS oc- 
cupation amplitudes {ipa, Va,a = 1, fi}. Thus the typ- 
ical mean-field state is a BCS state 
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where |0) is the vacuum state and Ua = y^l — v^. The 
product runs over all pairs of time-reversed partners 
(a, a) indicated by a > 0. The mean field equations 
are obtained by variation of the total energy with re- 
spect to single-particle wave functions and pairing am- 
plitudes. The expression for the total energy is the key 
ingredient in the modeling. Wc will here employ the 
Skyrme energy-density functional together with a pair- 
ing functional. Both functionals are expressed in terms 
of a few local densities and currents: local density p, 
kinetic-energy density t, spin-orbit density J, current j, 
spin density cr, kinetic spin density r, and pair current 
^. All appear twice, for protons and for neutrons, e.g.. 
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It is often useful to recouple to sum and difference, e.g., 

p^ Pp + Pn , p = Pp- Pn , (2h) 

and similarly for all other densities and currents. The 
sum plays a role in the isoscalar terms of the energy func- 
tional and we will call it the isoscalar density p. In a sim- 
ilar manner, the difference plays the role of an isovcctor 
density p. 

Our starting point is then the most general Skyrme 
energy functional 
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Accounting for the slight difference between the proton 
and the neutron mass (ji^ jlm^ = 20.749821, ?i^/2m„ = 
20.721260) becomes important for exotic and heavy nu- 
clei. The B {B') parameters determine the strength of 
the isoscalar (isovector) forces. The principle Skyrme 
functional £skyrmc contains just the minimum of time- 
odd currents and densities which is required for Galilean 
invariance [2^ . namely the combinations pr — and 
pV-J 4- cr-(Vxj). Further conceivable time-odd couplings 
play a role only for odd nuclei and magnetic (unnatural 
parity) modes [13 • These additional time-odd terms 
are often derived starting with a density-dependent zero- 
range force. We take here the point of view of density- 
functional theory and start from an energy functional, 
choosing the minimalistic approach in the time odd chan- 
nel. That modeling is similar to the strategy which is also 
pursued in the relativistic mean-field model [ll|, [l^, [13] . 
The parameters Bi , B'^ are more convenient for the func- 
tional form ([3]). They are uniquely related to the widely 
used standard Skyrme parameters ti,Xi through 
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Note the two non-standard entries. There is an extra 
parameter 64 for tuning the isovector dependency of the 
spin-orbit interaction. A zero-range two-body spin-orbit 
interaction leads to the fixed relation b'^ — t4/2 while 
a spin-orbit structure as suggested by the RMF is asso- 
ciated with 64 = [1^. We use it here as free parame- 
ter. The tensor spin-orbit term associated with the terms 
oc Ci,C[ is omitted in many parameterizations. Here we 
introduce a new parameter 77tis as a switch factor where 
Tytis = 1 includes the full tensor spin-orbit and Tytis = 
selects the widely used option to ignore the tensor spin- 
orbit term. The Coulomb functional (j3dp depends only 
on the charge density and stays outside this distinction. 
Its second term approximates exchange in the Slater ap- 
proximation [1^ . Note that we use the proton density in 
place of the charge density. This is a widely used, more or 
less standard, approximation. The center-of-mass correc- 
tion (|3f| is an approximation to the full cm. projection 
[2?! . It is applied a posteriori and its contribution to the 



mean-field equations is neglected. There are other recipes 
for the cm. correction not considered here, see [3, [iBl- 



The pairing functional ([3c|) involves the pair current £, 
and can be derived from a density-dependent zero-range 
force. The parameter ppair determines the weight of den- 
sity dependence. The limit ppair — > 00 recovers the 
pure S interaction (DI) which is also called volume pair- 
ing. The general case is the density-dependent 6 inter- 
action (DDDI). A typical value near matter equilibrium 
density ppair = 0.16 fm~'^ concentrates pairing to the sur- 
face. Thus it is often denoted as surface pairing. We will 
consider ppair as a free parameter and it will lead to an 
intermediate stage between volume and surface pairing. 

The results of the BCS calculations depend on the 
space of single-nucleon states taken into account, called 
here pairing phase space. In fact, the cut-off is 
part of the pairing description. It is provided by 
the phase-space weights Wa in the above pairing func- 
tionals. We use here a soft cutoff profile such as, 
Wa = [1 + exp((ea - (ef + ecut))/Ae)]"^ where typi- 
cally ecut = 5 MeV and Ae = ecut/10 [li,[li]. This works 
very well for all stable and moderately exotic nuclei. For 
better extrapolation ability away from the valley of sta- 
bility, the fixed margin Ccut may be modified to use a 
band of fixed particle number cx A^^/s instead of a fixed 
energy band [30| ■ 

When checking the performance of the fits for vari- 
ous observables (see section IIIIPp . we encounter also 
deformed configurations. These are computed with a 
code allowing for axially symmetric and reflection asym- 
metric configurations. The reflection asymmetry be- 
comes important for fission barriers in actinides. For 
well deformed nuclei, the dominant part of the cor- 
relation energy comes from angular-momentum pro- 
jection. We have accounted for that at the level 
of the Gaussian-overlap approximation (GOA) in a 
form which provides correctly a smooth transition to 
spherical shapes where the correction vanishes [U, 
[33 . [H, [ll]. This amounts to subtract i?zPE,rot = 
5(('^^)/4)('^^)/(20iot) where the interpolating function 
is g{x) ~ xdxlog[Ji^ daexp{x{a^ — 1))) . That cor- 
rection is included in all following results dealing with 
deformed configurations. 



III. FITTING STRATEGY 

A. Global quality measure and minimization 

The free parameters of the SHF ansatz arc going to 
be determined by a least squares fit. To that end, we 
build a global quality measure by summing the squared 
deviations from the data as 
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where O stands for one of the selected observables, "nucl" 
for a nucleus, the upper index "th" for a calculated value, 
and "cxp" for the experimental date. The denominator 
AO stands for the adopted error of that observable. It 
renders each contribution dimensionless and regulates the 
relative weights of the various terms. The experimental 
uncertainty is of little help here because the experimental 
precision of these basic bulk observables is much better 
than what we can expect from the mean-field description, 
particularly for the binding energy. The limiting factor 
comes from theoretical uncertainties, i.e., the quality we 
can expect from a mean-field description. We will briefly 
describe the observables in section lTlI Bl and the selection 
of nuclei together with the adopted errors in section lTIl CI 

The total quality measure is a function of all model 
parameters, i.e. — ■■■,Pn)- We search for those 

parameters which minimize x^- standard technique 
we employ the minimization technique from Beving- 
ton [3a, [3al , complemented occasionally by Monte-Carlo 
sampling to enhance the chances for ending up in the 
global minimum. 

The aim of this publication is to explore systematically 
the various influences from key features as, e.g., incom- 
pressibility or symmetry energy. We thus add constraints 
on such key features. That is done most simply by adding 
the wanted features as additional observables with very 
small adopted errors AO to the x^- 

The rules of x^ fitting also provide information to es- 
timate the statistical errors for extrapolations to other 
observables. Let us consider some observable A. Its ex- 
pectation value is a function of the model parameters, i.e. 
A = A{pi, ...,pn). The extrapolation error then becomes 
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That is the allowed variation of the observable A within 
the ellipsoid of x^ ~ Xmin — 1j i-*^- ^^r all x^ which stay 
at most one unit above the minimum. We will exploit 
that feature to compute extrapolation errors for all ob- 
servables not included in the fit data. 



B. Selection of fit observables 

The total binding energy Eb of a nucleus is the most 
immediate observable in self-consistent mean-field mod- 
els. It is naturally provided by the numerical solution 
of the mean-field equations. The next important observ- 
ables are related to the spatial extension of the nucleus. 
It can be accessed experimentally through elastic electron 
scattering which yields the nuclear charge form factor 
Fchiq) and by Fourier transformation the nuclear charge 
density Pch{r) H^]. SHF calculations yield first the pro- 
ton and neutron densities which need to be folded with 
the intrinsic electro-magnetic structure of the nucleons 
to obtain the charge density and charge form factor [s^ . 
Mean-field models provide reliable information about the 



form factor at low momentum q [H, . This region of 
the form factor at low momentum can be described by 
three parameters [st} . the r.m.s. radius 



i^ch(O) dq- 



ri^ch(g) 



q=0 



the (first) diffraction radius 



R = 



4.493 

(1) ' 
% 



(6) 



(7) 



which is determined from the first zero of the form factor 
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0, and the surface thickness 
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Therein i^box('?) corresponds to the form factor of a ho- 
mogeneous box of radius R, i.e. i^box('i') — '^jiiqR) / [qR)- 
The diffraction radius parameterizes the overall diffrac- 
tion pattern which resembles that of a box of radius R. 
It is the box- equivalent radius. For more details find a 
recent summary in 

The most prominent observable for T = 1 pairing cor- 
relations on the mean-field level is the odd-even stagger- 
ing of nuclear masses from which an approximation to 
the pairing gap can be extracted, e.g., using a five-point 
formula yields for the neutron gap 

A(,5) = -^E{Z,N + 2) + ]^E{Z,N + l)--^E{Z,N) 



^]^E{Z,N -l)-^-E{Z,N^2) 



and similarly for the proton gap A^p \ The computation 
of odd nuclei using mean-field models, however, is very 
involved. It requires deformation and a long search for 
the optimum blocked configuration which increases nu- 
merical expense by two orders of magnitude. Thus it is 
not well suited for systematic surveys. One can compute 
an average pairing gap from the state-dependent gap Aq, 
using a weight which is sensitive to the region about the 
Fermi surface HEll,!!! 



(9a) 
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These spectral gaps A^ are found to be fairly well related 

to the five-point gaps A^^^ in mid-shell regions [sol. [i^. 
For magic nuclei and next to them, these two quantities 
develop very differently [i^. There remains some influ- 
ence from nuclear shape fluctuations on Aq^^ [ij which 
gives an uncertainty of about 10-20% for the comparison 
of the two different definitions ^ for a gap. 

We also need some information which is specific to the 
spin-orbit terms in SHF. The obvious quantity for the 
purpose are the l*s splittings of single-particle energies. 
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Experimental information on single-particle energies of 
even-even nuclei is drawn from the single-nucleon removal 
energies or from the low-lying excited states of the ad- 
jacent odd-A nuclei. That identification requires that 
the polarization effects induced by the extra nucleon (or 
hole) are small. The magnitude of these effects have been 
investigated for doubly-magic nuclei by [45| within the 
RMF and [i^ with SHF using the linear response the- 
ory. It is found that energy differences amongst particle 
states and hole states separately are robust. Thus spin- 
orbit splittings are a robust signal as long as they do not 
cross the shell gap. Furthermore, the involved single- 
particle energies should not be too far away from the 
Fermi energy to keep perturbation from higher configura- 
tions (additional Iph couplings) low. Another limitation 
has to be considered for states in the particle spectrum 
because the state should stay safely off the particle con- 
tinuum. Thus we decide to include only spin-orbit split- 
tings of hole states close to the Fermi energy in doubly- 
magic nuclei. The chosen data are listed in tables [1111 and 



C. Correlation effects and adopted errors 

The goal is to adjust an effective nuclear energy-density 
functional. The observables and nuclei included in the 
adjustment should be well adapted for a pure mean- field 
description. Thus we have to scrutinize the most dom- 
inant effects going beyond mean field. As pointed out 
in the introduction, the strong short-range correlations 
as well as correlations from resonances and high modes 
follow a trend smooth in nucleon numbers Z and N and 
can be assumed to be effectively incorporated into the 
energy-density functional We have to care, however, 
about correlation effects which vary strongly within the 
chart of isotopes. These stem predominantly from the 
low-lying quadrupole excitations, associated with large 
amplitude collective motion like soft vibrations and rota- 
tions, for extensive recent analysis see [l^ [10, [U • We are 
concerned here with the influence of correlations on bulk 
observables, the change in energy SEcon- = -E-corr ~ E^i 
(where -Ecorr is the energy from the correlated calculation 
and i?nif the mean field result), and similarly the changes 
6r in charge r.m.s. radii, SR in diffraction radii, and 6a 
in surface thickness. 

A good chance for small correlation effects exists for 
semi-magic nuclei where either the proton or the neutron 
number corresponds to a shell closure. These nuclei arc 
generally spherical. Figure [T] summarizes the correlation 
effects on binding energies (lower panels) and radii as well 
as surface thicknesses (upper panels). The results were 
obtained with the parameterization SkI3 |25| . Other pa- 
rameterizations yield very similar results |2l| such that 
figure[T]is a typical result for the error distribution for any 
reasonable Skyrmc force. The correlation effects gener- 
ally shrink with increasing system size, which is expected 
because low-lying excitation energies generally decrease. 



It is a bit surprising, however, that the correlation en- 
ergies in isotopic chains can grow so large in the mid- 
shell region. The hitherto often underestimated isotonic 
chains are visibly less perturbed. Light nuclei generally 
show larger fluctuations, sometimes even acquiring un- 
physically positive correlation energies. The small posi- 
tive values about 0.2-0.5 MeV are still within the preci- 
sion of our method and imply practically negligible corre- 
lation energies. The only exception is ^°Ca at (Z=N=20). 
The unreasonably large value is due to the fact that the 
2+ mode is not really collective in that nucleus [2l[, but 
that we compute correlations using the GOA which is not 
necessarily valid for non-collective modes. This unphys- 
ical case, however, is anyway excluded from the fit data 
set due to the Wigner effect (for iV = Z) [irj . 

The faint horizontal lines in figure [1] indicate the error 
bands associated with the adopted errors, 1 MeV for Eb, 
0.04 fm for i?diff and cr, and 0.02 fm for Tims- Observ- 
ables which stay below these limits are included in the fit 
data. Points which are outside the desired error bands, 
but not too far away, are included with degraded error 
weight which is done by multiplying the general adopted 
error AO for that nucleus with a certain factor depending 
on how large the expected correlation could be. More- 
over, all nuclei with N=Z are excluded already because 
they carry a correlation contribution from the Wigner 
energy [47| and we do not yet have reliable means to 
compute that correction. The choice of fit data thus de- 
duced is listed in tables IIIII and IIVI One finds therein 
also the values for pairing gaps and spin-orbit splittings. 
Wc do not have any reliable estimate for the errors in 
these observables. The pairing gap is deduced from the 
five-point difference of binding energies [30j . The pairing 
gaps for nuclei whose energy has degraded error weight 
are modified by the same degradation factor for reasons 
of consistency. The overall error is tuned such that the 
r.m.s. average x from all pairing gaps is of the same or- 
der as for the other observables. A similar strategy taking 
over weight factors from binding energies is applied for 
the spin-orbit splittings. Altogether this amounts to an 
adopted error of 0.12 MeV for the gaps and 10% for the 
spin-orbit splittings, for details see again tables IIIII and 

IIYl 

Correlation energies are always negative (the few ex- 
ceptions at low A are defects) and r.m.s. radii always 
grow, but the fits produce an average with deviations ly- 
ing on both sides. One would very much like to leave 
a margin for correlations. A safe construction would re- 
quire a deeper knowledge of the systematics of correlation 
effects and development of simple estimates of them. At 
the present stage, we to fit to straightforward mean-field 
results and refrain from adding empirical corrections. 



D. More detailed observables 

There are many more nuclear properties which are in- 
teresting to look at but which arc not yet suited for in- 
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chains of semi-magic isotones chains of semi-magic isotopes 

Z=20 Z=28 Z=50 Z=82 N=20 N=28 N=50 N=82 N=126 



r.m.s. radii 5r 




20 28 28 4050 66 82 114 128 20 20 28 40 5050 66 82 92 
Neutrons N Protons Z 



FIG. 1; Correlation effects on semi-magic isotopic (left) and isotonic (right) chains computed with SkI3 and DI pairing. The 
isotopic chains are Z=20, Z=28, Z=50, and Z=82, and the isotonic chains N=20, N=28, N=50, N=82, and N=126. Upper 
panels: correlation shifts of charge r.m.s. radii (Sr), diffraction radii, (SR), and surface thicknesses {Sa). Lower panels: 
correlation energies (SE^ott)- Horizontal dotted lines indicate intended error limits. 



elusion in a systematic fit, beeause the relation to ex- 
perimental data is somewhat uncertain (neutron radii, 
giant resonance frequencies, fission barriers) or too cum- 
bersome to compute (low lying collective states). We will 
look at several such observables a posteriori: the neutron 
skin in ^°®Pb, the isotope shift of charge radii between 
^^*Pb and ^°*Pb, the neutron level sequence near the 
Fermi surface in ^^^Sn, extrapolation to super- heavy ele- 
ments (SHE), and giant resonances (GR). Some of these 
observables, namely the GR resonance peak frequencies, 
are even used as additional selection criteria. We will 
now explain their computation, while the other observ- 
ables are introduced in later sections. 

The dominant excitation modes of the nucleus are 
the giant resonances (GR). Their average peak posi- 
tion can be related to basic features. Heavy nuclei 
show least spectral fragmentation and are best suited 
for evaluating these averages. We will consider GR 
in ^°^Pb, in particular the isoscalar giant monopole 
resonance (GMR), the isovector giant dipole resonance 
(GDR), and the isoscalar giant quadrupole resonance 
(GQR). The spectral strength distribution is computed 
by the Random-Phasc-Approximation (RPA) done self- 
consistently with the same Skyrmc interaction as was 
used for the ground state, for technical details see 
[39I lisj . We use a large phase space on a large spher- 
ical grid of 30 fm radius to achieve a sufficiently fine 
discretization of the continuum [ioj for subsequent fold- 
ing with a Lorentzian of frequency dependent width 
F = max((;iw-8MeV)/3.5MeV,0.1MeV). The linear 
uj dependence starts at neutron emission threshold with 



an empirically adjusted slope. It simulates the escape 
width and to some extent the coUisional width. The 
strengths in ^°*Pb all have one unique peak in the GR 
region frequency that can easily be read off. Compar- 
ison with the experimental data will be done with re- 
spect to that GR peak frequency. The experimental val- 
ues are Kujgmr ~ 13.7 MeV, huJcpR ~ 13.6 MeV, and 
?iwGQR = 10.9 MeV [13, mi, Hi. 



E. Nuclear matter properties (NMP) 

Homogeneous nuclear matter describes the leading 
contributions (volume terms) to nuclear properties. They 
are considered as useful pseudo-observables characteriz- 
ing the bulk properties of effective interactions. There 
exist close relations between such bulk properties and cer- 
tain combinations of Skyrme parameters, but it is often 
more instructive to discuss a parameterization in terms 
of these nuclear matter properties (NMP). When per- 
forming systematic variations of forces, we will consider 
a scan of dedicated values for selected NMP rather than 
simple Skyrme parameters. 

The leading quantity is the energy per particle E/A{p), 
often called equation of state. Its minimum {E/A)cq at 
saturation density pcq defines the volume energy, related 
to the equilibrium state of nuclear matter. Energy and 
density in finite nuclei are modified by surface and shell 
effects but always stay close to these guiding values. Fur- 
thermore, we discuss NMP related to excitations (zero 
sound). The incompressibility at the saturation point is 
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given by 



9p' 



cq 



(10) 



where p = pn + Pp is the total density. It corresponds to 
the curvature at the minimum, and is related to breath- 
ing modes like the giant monopole resonance [s^. The 
symmetry energy coefficient is related to the isovector 
curvature at the saturation point 



1 



E 



-''sym 



2 d(p„ - pp)2 A 



(11) 



cq 



Finite nuclei arc also sensitive to smaller densities whose 
symmetry energy is characterized additionally by the 
slope 



E 



2 d{pn + Pp) d{pn - PpY a 



cq 



The isoscalar effective mass m* is calculated as 



2m* 



2m 



d_E 
d^A 



(12) 



(13) 



cq 



The isovector effective mass is usually expressed as the 
enhancement factor of the Thomas-Reiche-Kuhn sum 
rule l54i. which reads 



2m 



d 



E 



ktrk 



d{Tn - Tp) A 



(14) 



Note the subtle difference between total derivatives in 
eqs. (jlOllip and partial derivatives in eqs. (|13ll4p . The 
latter take E as written in the Skyrme functional ^ con- 
sidering all densities and currents (p, t, ...) as indepen- 
dent variables while the total derivative first expresses 
all quantities in terms of the actual Fermi momentum 
k-p = kp{p), and thus the density, before performing the 
derivative, e.g., considering r — > t{p). 



IV. RESULTS AND DISCUSSION 

In this section, we will present the parameterizations 
obtained by adjustment to the above selected set of 
data, optionally with an additional constraining condi- 
tion. The aim of the survey is to explore the influence 
of varying conditions in a systematic manner. We thus 
name the forces with the header "SV" and add differ- 
ent three digits qualifiers to indicate the constraint for 
which a force was adjusted. As outlined in section|TTl the 
Skyrme energy functional leaves some options to choose. 
The practical consequences thereof will be discussed in 
section IIVDI As a result, we will employ the following 
standard choices: The tensor spin-orbit is omitted, i.e. 
77is = 0, while isovector spin-orbit coupling 64 and the 
cutoff density in the DDDI recipe Ppair are allowed as 
free parameters in the fits. 



Furthermore, we will often compare with a few typical 
parameterizations from the existing literature: SkM* as 
a widely used, meanwhile somewhat obsolete, old stan- 
dard [Sal . It belongs to the second generation of Skyrme 
forces which for the first time delivered a high-precision 
description of nuclear ground states (as compared to the 
first generation forces). It was developed with an explicit 
study of surface energy and fission barriers in scmiclas- 
sical approximation. The set SLy6 and its cousins have 
been developed with a bias to neutron rich nuclei and 
neutron matter aiming at astrophysical applications [5^ . 
SkI3 (and SkM) exploit the freedom of an isovector spin- 
orbit force to obtain an improved description of isotopic 
shifts of r.m.s. radii in neutron rich Pb isotopes which 
posed a severe problem to all conventional Skyrme forces 
[25j . Finally, we consider one representative of the 
series of forces started in [l3|- These forces are fitted 
predominantly to binding energies, but employ an huge 
pool of nuclei including odd and deformed ones. For de- 
formed nuclei a simple correction for the angular momen- 
tum projection was added, and an ad-hoc correction for 
the Wigner energy in N^Z nuclei [IJ was applied. All 
these forces are are of comparable quality with respect to 
the reproduction of the binding energies of finite nuclei. 
They differ in details of treatment or boundary condi- 
tions. We choose here the force BSk4 from [s^l because 
it has an effective mass of m*/m = 0.92 which comes 
close to the typical values of our fits. 



A. Unconstrained fits 

1. The fit to standard data: SV-min 

In a first round, we follow the most straightforward 
strategy and adjust the model parameters to the data as 
selected in the previous section and detailed in tables Hill 
and IIVI without any additional constraint. The result 
is called the parameterization "SV-min" , for its detailed 
model parameters see table [V] Figure ^ shows the de- 
viation from the given data for each observable and nu- 
cleus. All results stay fairly well within the chosen error 
bands. That holds even for the points fitted with lower 
weight (open circles). In fact, the error bands are not 
fully exhausted and the r.m.s. errors stay safely below 
the adopted errors. The bands are generally filled on 
both sides of the zero line which means that the fit aver- 
ages nicely through the deviations. An exception is here 
the diffraction radius (second panels from below) where 
the deviation is always positive which indicates that this 
observable is not easy to adjust within the given model. 
The theoretical results are (within the allowed errors) 
systematically larger than the experimental data. Only 
few points exist for the spin-orbit splitting (uppermost 
panel) and these do not fit as nicely as the other ob- 
servables even with a low demand such as a 20% r.m.s. 
error. This indicates that single particle structure is a 
very demanding observable. We will see that again for 
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FIG. 2; Deviation from experimental data for the nuclei and observables in the fit sample as listed in tables ITITl and ITVl The 
results are drawn versus A'' for the isotopic chains (left) and versus Z for the isotonic chains (middle). The values (perfect 
matching) and the adopted error bands are indicated by dotted horizontal lines. Filled circles indicate data points fitted with 
full weight while open circles stand for data points with reduced weight (weight factor > 1 in tables Hill and HV)) . Results are 
shown for the parameterization SV-min resulting from a straightforward, unconstrained fit to the data. 
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x' 


/point 


r.m.s. error 


binding energy Eb 


12.07 


0.17 


0.62 MeV 


(1.0) 


diffr., radius R 


11.18 


0.40 


0.029 fm 


(0.04) 


surface thick, a 


4.22 


0.26 


0.022 fm 


(0.04) 


r.m.s. radius r 


15.86 


0.32 


0.014 fm 


(0.02) 


pairing gap Ap 


4.27 


0.25 


0.11 MeV 


(0.12) 


pairing gap A„ 


2.43 


0.15 


0.14 MeV 


(0.12) 


1-s splitting 


3.18 


0.45 


0.25 % 


(20) 


total 


53.22 


0.26 







TABLE I: Global quality measures for various classes of ob- 
servables as achieved with the parameterization SV-min. The 
second column shows the contribution from an observable 
to x' while the third column expresses this as P^r data 
point. The last column produces the r.m.s. errors as such 
and the numbers in brackets indicate the adopted error taken 
as weights for the fit, see eq. (|4]). 



other level sequences later on. 

The global quality measures for the fit SV-min are 
shown in table HI The r.m.s. errors (last column) show 
again what we have seen in figure[21 namely that energies 
and form parameters (radii, surface thickness) perform 
very well, better than the adopted errors as estimated 
from expected correlation effects. It seems that a part of 
the correlations can be accounted for by the model pa- 
rameters. Pairing gaps are just at the wanted limits and 
the 1-s splittings are somewhat at the edge. The good 
overall performance yields a very low total as can be 
seen from the very small per data point of about 1/4. 
Typical fits aim at a value of one. The strict rules of x^- 
fitting would allow to reduce the adopted errors until a 
X^ per data point of one is reached. We are not pursuing 
that strategy here because our adopted errors are deter- 
mined by the expected reliability of mean field models 
from estimating correlation effects. The fact that x^ per 
data point comes out much lower than one indicates the 
enormous versatility of the Skyrme energy functional to 
describe global ground-state properties. 

In order to check the interpolation and extrapolation 
properties, we show in figure [3] the errors in binding en- 
ergies and r.m.s. radii throughout all known nuclei in- 
cluding deformation. The energies are displayed in two 
ways. The lowest panel shows the energies as they result 
straightforwardly from the mean-field calculations (in- 
cluding cm. correction). The fit nuclei (filled squares) 
perform very well while large deviations can develop for 
other nuclei. Note that the majority of deviations is 
positive which indicates that additional binding through 
correlation energies would correct in the desired direc- 
tion. For well-deformed nuclei, the dominant part of the 
correlation energy comes from angular-momentum pro- 
jection. We have accounted for that at the level of the 
Gaussian-overlap approximation (GOA) [2l|, [3l|, [s^ • The 
results with angular-momentum projection are shown in 
the middle panel of figure [31 The deformed nuclei up to 
Pb now perform very well. The remaining discrepancies 




fit nuclei 



50 100 150 200 250 
mass number A 

FIG. 3; Errors in binding energies (lower two panels), r.m.s. 
radii (second from above), and proton deformation through- 
out all nuclei for which experimental data were available. The 
lowest panel uses the pure (deformed) mean-field value for 
the binding energies while angular-momentum projection is 
accounted for in the middle panel. The experimental defor- 
mations are deduced from electro-magnetic B(E2) values [s^ 
and are compared with the quadrupole variances (see section 
IIIID[) . For energies and radii, the nuclei which were included 
in the fit are marked by filled squares, well-deformed nuclei 
by open circles, and all others by triangles. 
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in the region A < 210 are very likely missing correlations. 
This means that the fit interpolates nicely for all these 
nuclei provided soft nuclei are computed with quadrupole 
ground-state correlations. However, for actinides and 
super- heavy elements (SHE), the trend of the deviations 
is too strong to be cured by remaining vibrational cor- 
relations. The extrapolation to deformed super-heavy 
elements is plagued by a growing trend to undcrbinding. 
We will take up that question later on. The results for 
the r.m.s. radii look agreeable. The larger deviations 
for some soft nuclei may still be cured by correlations. 
There are not yet enough data to read off a trend for 
super-heavy elements. The uppermost panel shows the 
difference between theoretical and experimental proton 
deformations as derived from the B(E2) values. They are 
given as dimcnsionless quadrupole moments associated 
with the operator (3 = Q20\/^ / (V^Ar^) where r is the 
r.m.s. radius and Q20 = r^Y2Q the spherical quadrupole 
operator. The theoretical values include the quadrupole 

variance, i.e. pf^^^'^ = = ^J0y + {A^$). It 

is to be remarked that this variance from the mean-field 
ground state underestimates the true variance for soft vi- 
brators and transitional nuclei. This lets us expect huge 
deviations for the majority of nuclei and this is indeed 
seen in the results. However, the picture looks better for 
rather rigid spherical nuclei with small correlations (our 
fit nuclei) and well deformed nuclei. It is a surprisingly 
nice agreement in view of the fact that sizeable contri- 
butions to the variance from the true collective ground 
state are still missing. 

Further properties of SV-min will discussed later in 
connection and comparison with other parameteriza- 
tions. 



2. A fit including a super-heavy nucleus 

We have seen above that the fitted force performs 
agreeably for interpolations, but not so well for extrap- 
olations. This suggests to extend the range of interpo- 
lation by adding super-heavy nuclei to the pool of data. 
Thus we have performed a fit where the binding energy 
of ^^^Hs was added with high weight to the standard set 
of data. That nucleus is well deformed. The angular- 
momentum projected energy has been taken as reference 
value for the fit. The results for the performance on en- 
ergies are shown in figure SI The error for ^^^Hs has suc- 
cessfully been curbed down to almost zero. But that is 
achieved at the price of lowering all energies for medium- 
heavy and heavy nuclei. The performance for the fit nu- 
clei is visibly degraded and other nuclei are turned to 
overbinding leaving no space for possible correlation en- 
ergies. We obviously encounter a deep-rooted problem 
with the given Skyrme energy functional when consider- 
ing the overall binding energy of SHE. The discussion of 
SHE will be continued in sections HVEel and HVETl 




50 100 150 200 250 
mass number A 



FIG. 4: Errors in binding energies (including cm. correc- 
tion and angular-momentum projection in GOA) throughout 
all nuclei for a parameterization where the binding energy of 
^®*Hs was additionally included in the x^- 



B. Fits with constraints on NMP 

1. Variation of NMP 

The uppermost entry of table Hi] shows the NMP and 
the final for SV-min together together with its ex- 
trapolation uncertainties computed according to cq. ([5]). 
The is very small in view of about 200 given data 
points. The NMP are more or less in commonly accepted 
ranges as used, e.g., in the liquid-drop model [s^ . The 
"ground-state" properties, equilibrium density peq and 
binding energy E/A, are well fixed while all the other 
NMP related to excitations show sizeable uncertainties. 
The fit leaves some freedom in these respects. Moreover, 
we will see that the most prevailing nuclear excitations, 
the giant resonances (GR), are not all so well tuned in 
SV-min. This suggests to exploit the freedom left by the 
fits for a fine tuning of GR. Moreover, it is interesting 
as such to explore the large space of still allowed varia- 
tions. To that end, we perform fits to the given set of data 
where additionally four NMP arc kept fixed: incompress- 
ibility K, effective mass 777,* /?7i, symmetry energy Osym; 
and sum rule enhancement factor n (see section fill El for 
its definition). The first two are isoscalar properties and 
the last two isovector. One may wonder why the slope of 
the symmetry energy, a'^y^, is not included in the varia- 
tion. The reason is that a'^y^ is tied up very closely to 
Osym by the fits such that only one of both shows that 
freedom of choice. 

A four-dimensional landscape of variations of NMP is 
too bulky to handle. We prefer to define one "base point" 
about which we perform variation of one NMP at a time 
thus dealing with four sets of variations. For the choice 
of the base point, we exploit the full space of variations 
in the four NMP and use the freedom to accommodate 
the GR properties as well as possible (see section ITV Cj) . 
After a longer search in that four-dimensional landscape. 
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force 


K m* /m asym k 


Pcq E/A Hsym 


2 

X 


SV-min 

± 


222 0.95 30.7 0.08 
8 0.15 1.4 0.40 


0.1610 -15.91 93 
0.0013 0.06 89 


53.2 


SV-bas 


234 0.9 30 0.4 


0.1596 -15.90 68 


57.9 


SV-K218 


218 0.9 30 0.4 


0.1615 -15.90 72 


57.0 


SV-K226 


226 0.9 30 0.4 


0.1605 -15.90 71 


56.3 


SV-K241 


241 0.9 30 0.4 


0.1588 -15.91 65 


61.3 


SV-maslO 


234 1.0 30 0.4 


0.1594 -15.91 59 


60.8 


SV-mas08 


234 0.8 30 0.4 


0.1597 -15.90 84 


59.5 


SV-mas07 


234 0.7 30 0.4 


0.1500 -15.89 109 


68.0 


SV-sym28 


234 0.9 28 0.4 


0.1595 -15.86 15 


63.4 


SV-sym32 


234 0.9 32 0.4 


0.1595 -15.94 119 


58.0 


SV-sym34 


234 0.9 34 0.4 


0.1592 -15.97 169 


61.6 


SV-kapOO 


234 0.9 30 0.0 


0.1598 -15.90 82 


52.6 


SV-kap20 


234 0.9 30 0.2 


0.1597 -15.90 74 


56.5 


SV-kap60 


234 0.9 30 0.6 


0.1595 -15.91 61 


60.5 


SV-tls 


234 0.9 30 0.4 


0.1595 -15.89 69 


61.2 


SkM* 
SLy6 
Skis 
BSk4 


217 0.79 30 0.53 
230 0.69 32 0.25 
258 0.58 35 0.25 
237 0.92 28 0.18 


0.1602 -15.75 95 
0.1590 -15.92 100 
0.1577 -15.96 212 
0.1575 -15.77 27 





TABLE II: NMP as defined in section lllTEl for the for the var- 
ious SHF parameterizations used in this paper {K in MeV, 
flsym in MeV, E/A in MeV, peq in fm"^, a'syn, in MeVfm^, 
m* /m dimensionless, and k dimensionless) . The rightmost 
column Usts the global quality measure x^- The parameter- 
ization SV-min results from an unconstrained minimization 
of the total quality measure according to eq. ([4]) with the 
data and adopted errors from tables |TIT] and |IV1 The other 
parameterizations were obtained by a fit constrained on four 
NMP. SV-bas is the base point of the systematic variation 
with the constraints; K = 234 MeV, m* /m — 0.9, asym = 30 
MeV, and k — 0.4. From that point, one property is varied, 
the incompressibility K (via power of density dependence a) 
in SV-K, the effective mass in SV-mas, the symmetry energy 
in SV-sym, and the sum rule enhancement in SV-kap. Fi- 
nally, SV-tls is constrained like SV-bas but employs the full 
tensor spin-orbit terms. Moreover, we append to the list the 
NMP for the four conventional Skyrme forces used in several 
comparisons. 



we find the following choice a good compromise: K = 234 
MeV, m*/m = 0.9, asym = 30 MeV, and k = 0.4. The 
corresponding parameterization is called "SV-bas" . We 
then vary each one of the four NMP while keeping the 
other three at the base value. This yields four sets of vari- 
ations with prefix SV-K for varied K, SV-mas for varied 
m*/m, SV-sym for varied agym, and SV-kap for varied 
K. Table HIl shows the NMP for these parameterizations 
and the detailed parameters for the functional ^ are 
provided in table |V] in appendix [X] There is also one 
parameterization SV-tls with the base values for NMP 
but now including tensor spin orbit, i.e 77tis = 1. This 
serves to explore the effect of the tensor spin-orbit term 



by comparison with SV-bas. 

The overall quality measure is shown in the right- 
most column of table [Tll The constraints, of course, 
degrade the quality a bit. But we see that the varia- 
tions yield which stay within an acceptable range of 
about 10% increase in x^- As we will see in section HV C[ 
SV-bas performs much better than SV-min for GR in 
208pj-,^ That counterweights the small losses on the side 
of ground-state properties. 

The NMP as given in table |TT] characterize ground 
states and excitation properties for modes with natural 
parity. As typical representatives for excitations with 
unnatural parity, we have also checked spin modes in nu- 
clear matter. These can be characterized by the Landau 
parameters go for pure spin excitations and for a spin- 
isospin mode [1^. Taking up the formula as given in 
[2^ , we have computed these Landau parameters for the 
energy functional Q . The only contribution comes from 
the tensor spin-orbit term. Thus we have go = = g'o for 
all forces with r/tis = 0. Only SV-tls has a non-vanishing 
tensor spin-orbit term. For that force we find go = —0.73 
and gQ = 0.191. These numbers stay safely above the 
critical value for spin(-isospin) instability. Thus all forces 
introduced here are stable in the spin channels. Note 
that we are here using the energy functional with min- 
imal time-odd terms, namely just those which are re- 
quired to achieve Galilean invariance. That leaves the 
effective interaction in all channels much more robust. 
Thus one should not be puzzled by the rather large neg- 
ative values for the spin-exchange parameters xi and X2 
in table [V] These come from expressing the functional 
in terms of the conventional Skyrme parameters, medi- 
ated through eq. ([3]). Spin stability is guaranteed in 
connection with the functional ^ A different picture 
would evolve when taking the Skyrme force literally as 
a zero-range force. This yields additional terms in the 
spin channel which can easily render a parameterization 
unstable in the spin channel. However, from an energy- 
density functional viewpoint, we see no compelling reason 
to include those terms. 



2. Trends of the errors 

The quality measure x^ is composed of different ob- 
servables whose relative weight is determined by the 
given adopted errors. It is known from earlier studies 
that the choice of observables has an influence on extrap- 
olated NMP (see, e.g., [1^). The sets with systematically 
varied NMP now allow to visualize these trends. Fig- 
ure [5] shows the r.m.s. averaged Xobs per data point for 
each observable separately. The figure has four panels to 
show the trends with respect to the four NMP variations 
considered: incompressibility K, effective mass m*/m, 
symmetry energy asym, and sum rule enhancement k. 
Variation of sum rule enhancement k changes very little 
in all observables. That feature is only loosely deter- 
mined by the data set, as already seen in table [Til from 
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FIG. 5: The distributions Xobs for each observable as indicated for four cases of systematically varied bulk parameters. Lower 
left: variation of symmetry energy asym; upper left: variation of sum rule enhancement k; lower right: variation of effective 
mass m*/m; upper right: variation of incompressibility. 



the rather large uncertainty for SV-min.. Wc will need 
further conditions to make a more definite choice. The 
other three features lead all to sizeable trends, but often 
lead in diiTerent directions. For example, the lower right 
panel shows that the r.m.s. radii would prefer low val- 
ues of m*/m while the energy prefers m* /m « 0.9, other 
observables to even higher m* jm. The final "optimum" 
for m* /m depends very much on the choice of the rela- 
tive weight of the different observables. Even within the 
energy as observable, we could revert the trend when giv- 
ing light nuclei more weight by using relative errors [36[ 
rather than absolute errors as done here. Significantly 
different trends are seen also for variation of agym and 
K. It is thus obvious that the relative weights in the 
composition of the determine the final extrapolated 
NMP. Note that the actual changes in the contributions 
to are small such that the total varies only very lit- 
tle when scrolling through the different NMP. There is a 
broad choice of well performing parameterizations. This 
explains why there are so many different SHF parameter- 
izations around which all provide a good description of 
nuclear ground-state properties but vary in several of the 
key features. In other words, the strategy of fitting 
leaves several vaguely determined aspects. One needs to 



include further observables which are more specific to the 
open features. 



C. Information from giant resonances 

1. GR in '^°^Pb and its relation to NMP 

Figure [S] collects results for the peak frequencies of 
the GR in ^ospb and of the GDR in ^^o. We concen- 
trate first on discussing the GR in ^"^Pb. The results 
for the straightforward fit SV-min are mixed. The GQR 
fits nicely, the GMR lies slightly to low, and the GDR is 
far off the goal. The uncertainties (see error bars for SV- 
min) are sufficiently large such that a good reproduction 
of GR in ^"'^Pb seems within the reach of allowed varia- 
tions. 

The relation between NMP and GR properties be- 
comes apparent from the various chains of the system- 
atically varied forces. The situation is particularly 
simple for the isoscalar excitations. There is a unique 
relation between an isoscalar GR and isoscalar NMP: 
the GMR is sensitive exclusively to K and the GQR to 
m* /m. Both these GR show a clean excitation spectrum 
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EGMR(^°^Pb) [MeV] Ec;DR(208pb) [MeV] EGQR('°^Pb) [MeV] ^gdrC^O) [MeV] 
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FIG. 6: Mean resonance frequencies for giant resonances in ^°*Pb (from left to right: GMR, GDR, and GQR) and for the GDR 
in ^^O computed with the various forces encompassing systematic variation of NMP properties incompressibility K, effective 
mass m* /m, symmetry energy Osym, and sum rule enhancement k as indicated at the right side of the panels. Additionally, 
results from SV-min, SV-tls and traditional Skyrme forces, as indicated, comprise the lower entries. Extrapolation errors for the 
three observables are shown as error bars for SV-bas and SV-min. Those for SV-bas apply to all forces fitted with constrained 
NMP. The experimental values are given in the uppermost entry and drawn through as faint vertical lines. 



with one peak. Thus we use these two data points to fix 
the otherwise weakly determined isoscalar NMP choos- 
ing m*/m = 0.9 to meet the GQR and K = 234 MeV 
to tune the GMR. The case is much more involved for 
the isovector GDR which reacts to two isovector NMP, 
to Osym and to K. Moreover, the spectral distribution 
(not shown here) tends to be strongly fragmented, par- 
ticularly for high m*/m and high k. Lower m*/m are 
excluded because we want to maintain the good adjust- 
ment of the GQR. Thus we stay with a compromise for 
the GDR, choosing k = 0.4 and a^ym = 30 MeV. This 
obviously does not perfectly meet the experimental peak 
position for the GDR. But the example of the GDR in 
^^O discussed later on demonstrates that the description 
of GDR by SHF is anyway not yet well under control. 
The compromise here is to be understood as a prelimi- 
nary setting, open for the necessary further studies on the 
GDR. Fixing these four settings in the fit yields the force 
SV-bas as introduced in the previous section. It serves as 
base point for further variations of NMP. The extrapola- 
tion uncertainties for SV-bas are, of course, much smaller 
than those for SV-min because the uncertainty in NMP 



has been fixed by choice. The strong reduction of the un- 
certainties confirms once more the close relation between 
GR and NMP. 

The force SV-tls is fitted as SV-bas, but with the ten- 
sor spin-orbit term switched on. There are only small 
changes as compared to SV-bas. These subtle shell ef- 
fects seem to have an only secondary influence on GR. 

The four lowest entries of figure [6] shows results from 
a few conventionally used Skyrme forces. The variation 
of the predictions is large confirming once more that GR 
are only loosely determined by ground-state fits and that 
explicit adjustment is needed for satisfying performance. 
Again, the GDR is not well described by any one of the 
four traditional forces. That does also hold for the force 
SGII which was developed originally for GR [6l|, and for 
which we obtain the GDR peak at 12.6 MeV, well within 
the results from other forces. 

2. GDR in the light nucleus ^^O 
The rightmost column in figure [5] collects results for 
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the peak position of the GDR in ^^O. Three of the four 
NMP show strong effects (k, m*/m, and agym)- The 
isovector chains along agym and k show the same trends 
which are at first glance natural in that the iJGDR(^^O) 
peak moves up with increasing i?GDR(^°^Pb). The sensi- 
tivity to variation of m* /m which is not present in ^"^Pb 
shows that the GDR in ^^O is more sensitive to shell 
effects than in ^"^Pb. But these are all comparatively 
moderate effects. The comparison with the average ex- 
perimental peak reveals a disaster. All results stay far 
below the goal. We have checked all conceivable varia- 
tions within the energy-density functional ([21), in earlier 
investigations (6^ and in the course of the present sur- 
vey, and did not find any way to come approximately 
close to the experimental i?GDR(^^0) without dramatic 
sacrifices on the quality of the ground-state description. 
Note that also none of the conventional Skyrme forces is 
able to reach approximately the wanted peak frequency. 
We conclude that there is no way to achieve a satisfying 
description of the GDR throughout all nuclei with the 
functional ([3]). There is an urgent need for a thorough 
investigation of that case. 



D. Fixing open options of the model 

The SHF functional ^ leaves a few options open con- 
cerning the spin-orbit model and the pairing functional. 
It is worth to check whether the measure can help 
deciding about preferred choices. To that end, we start 
from SV-bas and vary each one of these options sepa- 
rately. That variation proceeds similar to the variation 
of NMP. We fix the selected parameter at a wanted value 
and refit the other force parameters again by minimiza- 
tion of with additionally constraining the NMP to the 
values as used in SV-bas. We so to say produce vari- 
ants of SV-bas with one more parameter (either Ppair or 
?7tcnsor-is Or 64) fixed. 



1. Pairing model: DI versus DDDI 

Figure [7] shows the x^ for variation of the cutoff density 
Ppair in the DDDI functional (pe|) . The results are drawn 
versus p~^^^ because that provides a better scale. The 
limits are: p~^:^ — *■ leads back to DI while p~^- = 

0.16-1 fm^ = 6.25 fm^ is the typical DDDI value. The 
minimum x^ obviously lies between these two limits and 
the gain is considerable. We thus decide to use ppair as a 
free parameter of the pairing model. 



2. The tensor spin-orbit term 

The spin-orbit model leaves the option to include the 
tensor term oc J" stemming from the kinetic interaction. 




1 2 3 4 5 6 
Ppair"^ P""^] 

FIG. 7: Total quality measure a^s function of the inverse 
density-switching parameter p~^i^ in the variable DDDI func- 
tional. The DI limit relates to p~\^ = 0. The standard DDDI 
parameter Ppair = 0.16 fm"'^ is indicated by a vertical dotted 
line. 
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FIG. 8: Total quality measure a^s function of the tensor- 
spin-orbit factor 77tonsor-is. The value »7tonsor-is = (indicated 
by a vertical dotted line) is associated with ignoring tensor l*s 
coupling and ?7tcnsor-is ~ 1 with the full term as derived from 
a zero-range kinetic interaction. 



Many parameterizations ignore that term. A recent com- 
pilation explored the impact of tensor spin-orbit with 
flexible weight [6^ without finding clear signatures for 
optimum values. We have introduced a continuous switch 
factor 77tis to allow inclusion as tunable parameter. Figure 
[8] shows the result of a variation of Tytis. There is a steep 
increase in x^ from rjus = to 77tis = 1 and an almost flat 
landscape extending towards moderately negative values 
of 7/tis. The extremely shallow landscape will make 
fits extremely cumbersome because there is no drive to a 
clear minimum and negative values for the tensor switch 
factor seem a bit unorthodox. We thus decide to freeze 
the option 7/tis = 0, i.e. omitting tensor spin-orbit, and 
to consider Tytis = 1 (= full tensor spin-orbit) occasionally 
as a separate option. 
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FIG. 9: Dependence of the total quality measure on the 
ratio between isovector l*s coupling 64 and isoscalar l*s cou- 
pling 64 = 14/2. 



3. Isovector spin-orbit term 

The other open point concerns isovector spin-orbit cou- 
pling. That is described by the parameter 64 in the 
Skyrme functional ([3|). A fixed isovector fraction with 
64 = t4/2 is the standard spin-orbit model in conven- 
tional SHF functionals. A non-relativistic limit from the 
RMF suggests 64 « as the appropriate choice [1^, [U ■ 
Figure [9] shows the function of the ratio 2b\/ti. 

A value of zero corresponds to the RMF preference and 
one to the standard spin-orbit model of SHF. There is 
a considerable sensitivity. The minimum comes close to 
the RMF situation for the present test cases. However, 
the optimal 64/64 depends somewhat on the NMP con- 
straints (e.g.. effective mass, sum rule enhancement), as 
seen in table [V] We thus decide to consider 64 as a freely 
fitted parameter of the model. 



E. Performance for other observables 

In the following, we will explore the effect of variation 
of NMP on various detailed observables. It often happens 
that only one NMP shows significant effects on a given 
observable. In such a case, we will show only the one 
most relevant variation. 



1. Neutron skin in Pb 

A conceptually simple observable is the neutron ra- 
dius. It complements the radius information gained from 
the charge form factor. Unfortunately, its experimental 
determination is model dependent because the strong in- 
teraction is involved [1^. With that precaution, we con- 
sider the neutron radius in ^°®Pb from 
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FIG. 10: The neutron skin, 



in ■^'^^Pb as computed 



it as neutron skin, i.e. as difference between neutron and 
proton radius, and take as experimental reference value 



by the parameterizations with varied symmetry energy flsym. 
The experimental value from [g^ is drawn as a faint horizontal 
line together with assumed errors of ±0.05 fm. The set SV- 
min is shown together with the extrapolation errors. The 
extrapolation errors for SV-bas are smaller than the symbol. 



^ 'p ~ 0.15 fm. The neutron skin turns out to depend 
exclusively on the symmetry energy agym while all other 
NMP have no effect at all. The trend is shown in figure 
1101 One may even add a result from any other force into 
that plot and all would line up nicely on the given slope 
[igI . (6^ . The minimal fit SV-min agrees with the data 
within the extrapolation error. On the other hand, one 
has to take into account that the experimental neutron 
radius is extracted employing model analysis of data [6^ . 
Some uncertainty is associated with that result. Wc have 
adopted an experimental uncertainty of 0.05 fm, which 
still leaves a huge degree of freedom for agym- Only the 
conventional RMF parameterizations with their rather 
large agym > 34 MeV can safely be excluded [3. It is 
highly desirable to have more reliable data for the neu- 
tron skin in heavy, neutron-rich, nuclei. That would pro- 
vide direct access to the symmetry energy. 



2. The isotope shift from ^°^Pb to ^^'^Pb 

Another observable related to nuclear shape is pro- 
vided by the isotope shifts of char ge r adii. These are 
well accessible by optical methods [63 and differences 
of radii put weight on aspects complementing informa- 
tion from radii as such. Particularly interesting here 
is the shift r'^C^^'^Ph) - r^C^^^Ph) which is usually not 
well described in SHF with the traditional form for the 
spin-orbit term but very well in RMF [2^, i68j] . We will 
look at that observable with respect to the experimen- 
tal result r2(2i4pb) _ r^(^"^Ph) = 0.6085 fm^ Fig- 
ure [11] collects results for that isotope shift. From the 
four NMP, m* /m shows the strongest effect, although 
the other NMP are not totally ignorable, the forces with 
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FIG. 11: The isotopic shift, r^(^^-^Ph) - r^(^°^Ph), of charge FIG. 12: The energy difference between two occupied neutron 
radii as computed by the parameterizations with varied ef- 
fective mass m*/m (full squares) . The experimental value is 
drawn as faint horizontal line [63] ■ Results from some stan- 
dard parameterizations are drawn as open squares and labeled 
with their names in the literature. The parameterization SV- 
tls is a new fit with constraints as SV-bas, but using the full 
tensor spin-orbit, i.e. ijtis ~ 1, see the SHF functional ([HJ. 
Extrapolation errors are indicated by arrows for SV-bas and 
SV-min. 



the chain of varying effective mass and some other Skyrme 
parameterizations as indicated. The experimental value is 
indicated by a faint horizontal line [70l | . Extrapolation errors 
are indicated by arrows for SV-bas and SV-min. 



varying K yield a variation of 0.02 fm^, agym spans 0.05 
fm^, and k 0.03 fm^. All these variations stay a fac- 
tor of two below the strongest one for m* /m. The non- 
negligible extrapolation error of SV-bas indicates that 
there are also other ingredients of the force which in- 
fluence the isotopic shift. The isovcctor spin-orbit force 
does so by construction and the density-dependence 
of the pairing force is supposed to also have some effect 
[69| . The difference between SV-tls and SV-bas indicates 
the large impact of tensor spin-orbit. The present study 
using systematic variation of NMP, however, reveals that 
the effective mass still has the leading influence. It was 
the low m* /m of the older generation of forces which en- 
hanced the impact of the isovector spin-orbit term with 
&4 such that it lead to the success of SkI3 and SkM in ad- 
justing the isotopic shift [1^. The present preference of 
high effective masses leaves little chance to come so close 
to the data as can be seen from the error bars on SV-min. 
This calls for new investigations on that subject. 



3. Neutron level sequence in ^'^^Sn 

The single-particle energies are a subtle observable in 
connection with mean-field models. Density functional 
theory does not give a guarantee that they are correctly 
described [l3| and indeed they are hampered by the self- 
interaction error. That error, however, leaves the energy 
differences untouched [Tlj and it is anyway less dramatic 
in nuclei. The robustness of single-particle energy dif- 
ferences has also been confirmed in the nuclear context 



[45l . |46[. Moreover, the detailed level structure plays a 
crucial role for the properties of SHE It is thus 

important to investigate the performance of mean-field 
models in that respect. Simple test cases can only be 
doubly magic nuclei to avoid perturbations from pairing 
and polarization effects. The levels are usually described 
fairly well in ^°^Pb to the extent that at least the level 
ordering is usually correctly reproduced and often also 
the detailed energy differences [J|. However, a particu- 
larly obnoxious case is the sequence of occupied neutron 
levels in ^^^Sn. The experimental ground state of ^^^Sn 
has a spin 3/2+ while all mean field models predict ll/2~ 
spin. We take the energy difference between these two 
states, £n,ihii/2 ~ £n,2d3/2i as One number characterizing 
the sequence. Figure [Tj] shows the results for the chain 
with varied m*/m. It is not surprising that this varia- 
tion has most effect because m* /m is closely related to 
shell structure. The next sensitive NMP is k (not shown 
here) which is related to the isovector effective mass. The 
other parameterizations in the figure gather between 1-2 
MeV and all theoretical results (including those from the 
RMF I, El) are far off the experimental value. 

The origin of this discrepancy is not clear. It may be 
related to a peculiarity of the high angular momentum 
state involved in that difference, either that its spin-orbit 
splitting is underestimated or that the mean position of 
high spatial angular momenta is too high. It would re- 
quire extensions of the SHF functional to cure one of 
these two features. But it could also be a problem with 
the interpretation of the excitation spectrum in ^^^Sn as 
neutron levels in ^"^^Sn. What has not yet been checked 
so far is the influence of particle-core coupling in ^'^^Sn. 
The problem of the level sequence calls for further thor- 
ough investigations. 
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FIG. 13: Performance for extrapolations to extremely 
neutron-rich Sn isotopes for the two parameterizations SV- 
min and SV-bas. Lower: Predicted extrapolation uncertain- 
ties for binding energies. Upper: Predicted extrapolation un- 
certainties for two- neutron separation energies S2n- 



A similar discrepancy is seen for the energy difference 
between the proton l/ill/2 and 2g7/2 levels in Sn iso- 
topes [zl] where all mean field models yet fail to repro- 
duce the isotopic trend of the splitting. Two states with 
high angular momentum are involved. This suggests that 
an insufficient description of spin-orbit splitting is the 
more likely source of trouble. But again, this has yet to 
be investigated in more detail. 



4-. Extrapolations to neutron-rich Sn Isotopes 



236 



u 



barriers [MeV] 
1 2 3 4 5 6 7 8 





isomer 


2. barrier 


1 . barr. 


exp 


■ 


ae 


■ 


SV-K21 


w 


p s 


E 


SV-K22 


m 


4> 


S 


SV-bas 


i 


O J 


; 


SV-K24 


■ 


i> ) 




SV-m07 




p 




SV-m08 


k 




/ ■ 


SV-bas 


'w 


^ } 




SV-m10 


m 


U 


. 


SV-S28 


p 


p 




SV-bas 


p' 


f 


; 


SV-S32 


p' 


f 


■ 


SV-S34 




p' * 


■ 


SV-l<00 


f 


p 


C 


SV-I<02 


• 


(b s 


; 


SV-bas 




5 


(~ 


SV-I<06 


p 


6 5 




SV-min 














SV-tIs 


p 


o * 




SItlVI* 


p 


o 


m ■ 


SLy6 


p 


o 




SI<I3 




o 




BSI<4 


p 


o * 





K 



m /m 



^sym 



FIG. 14: Axially symmetric fission barriers (star — first bar- 
rier, open circle — second barrier) and position of the isomeric 
minimum (filled square) in ^^^U for the variety of parameter- 
izations as indicated. The first barrier Extrapolation errors 
are shown for SV-bas and SV-min. The errors for SV-bas 
are representative for all forces with constrained NMP. The 
experimental values are taken from [t^. [TSj. 



Astrophysical applications for r-process nuclei involve 
extrapolations deep into the regime of neutron-rich iso- 
topes. Figure [T3] shows the extrapolation errors as esti- 
mated from the least-squares techniques for binding ener- 
gies and two-neutron separation energies along the chain 
of exotic Sn isotopes. The lower panel for binding ener- 
gies shows a systematic growth of the uncertainty when 
moving away from the valley of stability. Note that the 
freezing of crucial NMP by tuning of giant resonances in 
SV-bas reduces the uncertainty by a factor of two. The 
upper panel for two-neutron separation energies basically 
behaves similarly, but the errors are much smaller than 
for the energies as such. Differences of energies probe the 
response properties and these are obviously a bit more 
robust. Note also the peak at neutron number iV = 92. 



There is a small sub-shell closure which is particularly 
sensitive to shell structure. The increase of uncertainty 
here indicates that an aspect of shell structure comes 
into play which is not so well determined as the general 
trends. 



5. Fission barriers in U 

Fission of actinides is a crucial nuclear process. It is 
usually characterized in terms of fission barriers. We de- 
duce the barriers from quadrupole-constrained SHF cal- 
culations considering axially symmetric shapes and al- 
lowing for reflection asymmetry. A series of constrained 
calculations for a broad range of deformations spans the 
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fission path and the associated deformation energy sur- 
face. All states along the path are well deformed. We 
thus compute the energy with the correction from angu- 
lar momentum projection in GOA, as explained in section 
im That energy curve for actinides displays the typical 
double-humped barrier fl^ with one isomeric minimum. 
From that we read off two barrier heights and the iso- 
meric energy. It is to be noted that the first barrier 
tends to go through triaxial shapes, not accounted for 
here. A further lowering by 0.5-2 MeV can be expected 
from triaxiality [zl, [T^I . 

Figure [TJ] shows results for fission isomer and fission 
barrier for ^^ey typical example for an actinide. 

The values fit generally well to the experimental data, 
particularly when considering some triaxial lowering for 
the first barrier. The strongest infiuence from NMP on 
the first barrier comes from the symmetry energy agym, 
some effect is added from m*/m, and K as well as k 
remain basically inert. The order is reversed for the sec- 
ond barrier where m* /m has largest impact and Ogym is 
secondary. These two NMP are the crucial handles on 
fission properties including the isomer. This reflects that 
fission emerges from a subtle interplay of bulk properties 
(here flsym) and shell effects (here m*/m). The impact 
of shell effects through the spin-orbit force is seen also in 
the step from SV-bas to SV-tls which yields a small, but 
non-negligible, lowering of the barriers. 

It is interesting to note that the variation of results for 
the conventional forces is larger than the variation within 
all SV forces. That indicates a large sensitivity of fission 
properties to fitting strategies. 



6. A known SHE: "^^^ Hs 

Some SHE have already been produced such that data 
are available for probing the predictive power directly. 
We consider here ^^^Hs as sample. The nucleus has a 
well deformed ground state. In fact, it belongs to an is- 
land of deformed shell closures [ll, [s^l. We discussed 
it, amongst others, already in connection with figures 
[3] and m Here, we continue with more detailed varia- 
tions. Figure [15] shows the binding energy of ^^^Hs for 
all forces in our sample. The test case embraces extrap- 
olation to SHE together with deformation effects. It is 
thus very sensitive and we find that all NMP have some 
influence, the strongest coming from m* /m. But none of 
the allowed variations of NMP brings the result in any 
way close to the experimental value. The same holds 
true for the conventional Skyrme forces. That feature 
was already observed in [tI] and it was also found that 
relativistic mean field models behave quite differently. 
Thus we have added a few results from relativistic mod- 
els, NL-Z NL3 [13, and PC-Fl These tend 
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to overbinding where SHE generally yields underbinding. 
The result indicates a deep rooted structural difference 
between these two classes of models, and possibly missing 
terms in both. 
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FIG. 15: The binding energy of ■^®*Hs for the variety of pa- 
rameterizations as indicated. The lowest three entries show 
results from relativistic mean-field models [tI, [t^ . Extrapo- 
lation errors are shown for SV-bas and SV-min. The errors 
for SV-bas are representative for all forces with constrained 
NMP. 



Reactions in the landscape of SHE are not so much 
determined by binding energies as such. Differences of 
binding energies are more crucial, particularly the Qa 
value which characterizes a-decay. Figure [H] shows Qa 
for ^^^Hs in comparison to the experimental value. All 
forces provide rather nice agreement with the experimen- 
tal value. There is general trend to about 0.2 MeV under- 
estimation. This may be due to the fact that the daugh- 
ter nucleus Z =■ 106 is softer than the parent Z — 108 
such that some additional correlation effect may come 
into play. Nevertheless, the agreement is much better 
than what one could have expected from the variations 
in the binding energies (see figure I15p . That indicates 
that energy differences can be predicted more safely. 
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FIG. 16: The Qa value of Hs for the variety of parame- 
terizations as indicated. Extrapolation errors are shown for 
SV-bas. Those for SV-min exceed the bounds of the plot. 



7. At the upper end a spherical SHE: Z,N = (120, 182) 

Further up, there are very interesting, yet unmeasured, 
SHE. In that regime, we consider as test case element 
Z,N = (120,182). at the upper edge of experimental 
feasibility in the regime of the spherical valley of stability. 
It resides on an upward extension of an a-decay chain 
recently detected [s^ [13, . This is a spherical nucleus 
due to the proton shell closure at Z = 120 and a a low 
level density for the neutrons at their Fermi energy (sgI . 
IstI . IssI ] . We consider for that element the binding energy 
Eb as such, the Qa value 

QaiZ,N) ^ EB{Z-2,N-2) + EBi2,2)-EBiZ,N) (15) 

which characterizes a-decay, and the fission barrier 
Bfis- The latter quantity is deduced from computing 
the potential-energy surface along the axially symmetric 
quadrupolc deformation path, for details see [89| . includ- 



ing again the correction from angular momentum projec- 
tion. Most SHE have only a single fission barrier. 

Figure[T7]shows the results for the hypothetical nucleus 
(Z/N)=(120/182). We concentrate first on the binding 
energy (left column). The span of predictions from 
conventional parameterizations is huge. The force SkM* 
is far off all other results which could be explained by 
two reasons: First, it is the oldest force in the sample 
and it was adjusted on a smaller data base available at 
that time 



55| 



Second, and probably more important, 
SkM* is the only force in the sample which uses a differ- 
ent recipe for the cm. correction, namely to take only 
the diagonal elements of the Pcm in eq. ((3f| which has 
dramatic consequences for the extrapolation to SHE [gO] . 
The more recent forces are grouped somewhat better to- 
gether, and the whole set of SV forces shows compara- 
tively little variation. The extrapolation error of SV-min 
is still smaller than the difference to conventional forces 
and SV-bas comes down to an uncertainty of ±0.8 McV. 
The difference in the extrapolation errors between SV- 
min and SV-bas shows the influence of the loosely fixed 
NMP. A better determination of NMP through GR does 
also improve the predictive value in the regime of exotic 
nuclei. 

The right column of figure [T71 shows the fission barriers 
for the SHE 120/182. These are sensitive to all details 
of a parameterization due to the subtle interplay of bulk 
properties (surface tension. Coulomb pressure) and shell 
effects. All NMP yield 1-2 MeV variation of the result. 
The sizeable extrapolation errors of about 0.8 McV for 
SV-bas show that there are other effects also at work. 
The fission barrier is determined by bulk properties (i.e. 
NMP) as well as shell effects [s^. The latter can be seen, 
e.g., by the large effect of the tensor spin-orbit force, see 
the difference between SV-bas and SV-tls. In spite of 
all these sensitivities, the more conservative extrapola- 
tion error of SV-min indicates that the predictions can 
be taken within about 1.2 MeV reliability. It is comfort- 
ing to see that most of the conventional parameteriza- 
tions stay also within these bounds. Thus we see that 
SHF notoriously predicts the fission barrier for 120/182 
around 6 MeV with angular momentum correction and 
7-8 MeV without. That remains in great contrast to the 
RMF where all predictions come out much lower (89| . 

There is a competition between a-decay and fission in 
the decay channels of SHE. The middle column of figure 
[17]shows the Q„ value for the SHE 120/182. This quan- 
tity is directly deduced from a difference of binding ener- 
gies. One would expect similar trends as for the binding 
energies (left column), but differences can suppress one 
trend and amplify another. That is what happens for 
the Qa- The trend with m*/m looses importance and 
the incompressibility K acquires more weight, but the 
variations are generally much smaller than for the fission 
barriers. Except for SkM*, all forces agree within ±0.4 
MeV. The fact that the extrapolation error for SV-bas is 
much smaller than for SV-min indicates that NMP are 
the leading dcterminators, in contrast to fission barriers 
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FIG. 17: Results of the various forces for the super-heavy element (Z/N) — (120/182). Left column: total binding energy. 
Middle column: Qa value. Right column: fission barrier. Extrapolation errors are shown for SV-bas and SV-min. The errors 
for SV-bas are representative for all forces with constrained NMP. 



where shell effects have mueh larger influence. 

V. CONCLUSION 

In this paper, we have performed a survey of the phe- 
nomenological adjustment of the free parameters of the 
Skyrme-Hartree-Fock (SHF) energy functional. The in- 
put data for the fits are taken from basic nuclear ground- 
state properties: energies, charge radii, surface thickness, 
selected odd-even staggering of energies, and some spin- 
orbit splittings. These data and appropriate nuclei are 
selected carefully such that they have small correlation 
effects. Thereby we consider only the strongly fluctu- 
ating correlations from low-lying collective quadrupole 
states. The smoothly varying correlations from higher 
excitations (two-body collisions and resonances) are sup- 
posed to be incorporated effectively in the SHF func- 



tional. The investigation of correlation effects led to a se- 
lection of semi- magic fit nuclei with large extension along 
isotonic chains but surprisingly short isotopic chains. All 
selected nuclei are spherical to avoid ambiguities from 
the handling of deformed minima (angular momentum 
projection). A quality measure is built from summing 
the squared deviations from the data with appropriate 
adopted error weights. The parameters of the SHF func- 
tional are optimized by least-squares fits with respect to 
the selected sample of data. The emerging r.m.s. errors 
in the fit observables stay well (by a factor two) below the 
adopted input errors from correlation effects which shows 
that the SHF model has some flexibility to incorporate 
part of correlations from the low-lying excitations. 

It is found that the minimization leaves some free- 
dom in nuclear bulk properties. In order to explore the 
space of well-fitted forces thus becoming available, we 
have performed a systematic variation of bulk proper- 
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ties, characterized in terms of nuclear matter proper- 
ties (NMP). This is achieved by adding to the selected 
ground-state data constraints on four NMP: incompress- 
ibility if, effective mass jTi^ jvn^ symmetry energy Ugyni, 
and sum rule enhancement factor k. These NMP are 
varied systematically to produce a set of forces with dif- 
ferent properties, all having about the same high quality 
concerning the nuclear ground-state properties in the fit 
data. The set of parameterizations thus obtained was 
used for a thorough investigation of the predictive power 
of the SHF model by looking at the results for several 
detailed observables in stable and super-heavy nuclei. 

We have used the quality measure to check some 
open options of the present ansatz for the SHF functional. 
There used to be the decision between volume pairing 
which corresponds to a simple zero-range pairing force 
and surface pairing which augments that with a density 
dependence such that pairing is basically switched off 
for densities near bulk equilibrium. We have allowed a 
flexible switching with a density parameter between bulk 
equilibrium (surface pairing) and infinity (volume pair- 
ing). It turns out that the optimum choice is just between 
these extremes. The SHF functional allows a free choice 
for the isovector spin-orbit force. A free variation is found 
to be advantageous as compared to the standard choice 
which is linked to the isoscalar value or to the choice 
of zero isovector term deduced from the RMF. There is, 
furthermore, the choice to include the tensor spin-orbit 
term which is related to the kinetic zero-range two-body 
interaction. The quality measure prefers a model with- 
out tensor spin-orbit. All three decisions, however, are 
related to changes in the total by about 5%. 

We have looked at the separate contributions of an ob- 
servable (e.g., energy) to the total ^ . The trends of these 
contributions with systematically varied NMP pull often 
in different directions. This means that the choice of rel- 
ative weights of one observable with respect to another 
decides on the final result. The trends are weak and thus 
small changes in weight can cause large drifts in final 
NMP of a least-squares minimum. This is a quite un- 
desirable element of arbitrariness in the adjustment pro- 
cedure. On the other hand, the weakness of the trends 
means that there are several features not so sharply de- 
termined from ground-state fits. We exploit the freedom 
to additionally adjust the peak positions of three decisive 
giant resonances in ^°*Pb: the isoscalar giant monopole 
resonance (GMR), the isoscalar giant quadrupole reso- 
nance (GQR), and the isovector giant dipole resonance 
(GDR). The GMR fixes the incompressibility and the 
GQR the effective mass. Both resonances can be well de- 
scribed by the model. The GDR is less conclusive in two 
respects. First, it is sensitive to two NMP, the symmetry 
energy and the sum rule enhancement factor. And sec- 
ond, it does not yet allow fully satisfying tuning due to 
the strong fragmentation pattern. Moreover, we find that 
the GDR in ^^O is far off the experimental peak position 
for all reasonable SHF parameterizations. The present 
forms of the SHF functional do not yet allow a proper 



description of the GDR throughout all nuclei. This case 
of the GDR still requires further investigation. 

In summary, we have built one force SV-min with min- 
imal from an unrestricted fit, another force SV-bas 
from NMP constrained fits with GR fine-tuning, a series 
of forces with systematically varied NMP, and finally a 
variant of SV-bas with tensor spin-orbit term. Using this 
set of trial forces, we have investigated the performance 
for other nuclei and the predictions for a variety of more 
detailed observables. The interpolation to nuclei within 
the fitted mass range, i.e. A < 220, perform fairly well 
when including angular momentum projection. The re- 
maining undcrbinding of 1-2 MeV for soft nuclei is re- 
solved by vibrational correlations. The extrapolation to 
heavier, deformed nuclei, however, develops large under- 
binding up to 4 MeV. The attempt to cure that by in- 
cluding super-heavy elements in the fit did not work out 
well. It sacrifices too much of the quality reached for sta- 
ble nuclei. We encounter here most probably a problem 
with an insufficient form of the present energy functional. 
The binding energy of spherical super-heavy nuclei is a 
rather robust quantity depending mainly on the effec- 
tive mass and other features influencing shell structure. 
Deformed super-heavy nuclei are more sensitive because 
of the deformation energy. The fission barriers in super- 
heavy elements is more sensitive depending on almost any 
feature of the force, all NMP and more detailed entries as, 
e.g., the tensor spin-orbit term. Nonetheless, all predic- 
tions lie within a band of ±1.2 MeV out of about 10 MeV 
which is a very comforting result for such a subtle observ- 
able. Even more robust are the Qa energies determining 
the rate of a decay, which are most strongly influenced 
by the incompressibility. This shows that differences of 
energies can behave differently from the energies as such 
(where the effective mass was most influential). The neu- 
tron skin in ■^''^Pb depends exclusively on the symmetry 
energy. It would be the ideal means to determine this 
crucial isovector property, but the measurement needs 
to be reliable and precise. An uncertainty of 0.02 fm 
in the neutron radius translates to an uncertainty of 1 
MeV in the symmetry energy. Practically all existing 
SHF parameterizations comply with the presently avail- 
able data within their large uncertainties. The isotope 
shift of charge radii between ^^^Pb and ^°®Pb depends 
sensitively on the effective mass and other features in- 
fluencing the shell structure. The recent fits all tend to 
rather large effective mass which, in turn, raises problems 
with accommodating that isotopic shift, even when al- 
lowing full freedom in the isovector spin-orbit term. The 
case, which seemed to be solved for previous parameteri- 
zations with low effective mass, is open again and requires 
further investigations. 

Altogether, we find that the SHF model provides an 
excellent description of nuclear bulk properties. Ground- 
state properties alone, however, leave several features 
weakly determined. Additional information from excita- 
tion or more detailed ground-state observables is required 
to fix all aspects of the SHF model. At the same time. 
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the added observables reveal some insufficiencies of the 
model which call for further investigations. Most urgent 
seems to be an extension of the spin-orbit model and a 
better description of dipole giant resonances. 
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TABLE III: Experimental data for the fits, part I: along isotopic chains. Each column stands for an observable as indicated. 
The second line shows the globally adopted error for each observable. That error is multiplied for each observable by a further 
integer weight factor which is given in the column next to the data value. 
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TABLE IV: Experimental data for the fits, part II: along isotonic chains. Doubly magic imclei which would fit both sequences 
are not repeated here. For further explanations see table Hill 
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TABLE V: The parameters of the energy-density functional (|3]) for the various SHF parameterizations used in this paper. All 
parameterizations use a proton mass of /2mp = 20.749821 and a neutron mass of ?i^/2m„ — 20.721260. For an explanation 
of the labels see table |II] 
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